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We study the effect of magnetic impurities in the d a ,2_ y 2-wave superconducting (SC) state of 
the two dimensional t — J model. The spin-orbit and the spin-exchange interactions are examined 
by treating the impurity as a classical spin. The Bogoliubov de Gennes equation derived within 
a slave-boson mean-field approximation is solved numerically at T = 0. The spin-exchange 
scattering induces spin-triplet p-wave SC order parameters near the impurity, while a SC state 
with broken time-reversal symmetry and a spontaneous current appears in the presence of the 
spin-orbit interaction. When both interactions coexist, it turns out that a state which carries a 
spontaneous spin current occurs. 
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The symmetry of the superconducting (SC) states in 
high-T c cuprates has been studied intensively, and now it 
is established that the state has a predominantly d x 2_ y 2- 
wave character with a possible mixture of an s-wave com- 
ponent due to the orthorhombic lattice distortion in some 
systerns.cP Recently the effect of nonmagnetic (Zn) and 
Magnetic (Ni) impurities in high-T c superconductors are 
of particular interest both theoretically and experimen- 
tally. This is because the effects of randomness in such 
unconventional superconductors can be quite different 
from those in conventional superconductors. 

The low-energy electronic states of high-Tc cuprates 
can he. described by the t — J model on a square lat- 
tice.BEP Mean-field (MF) theories based on a slave-boson 
method predict a superconducting state with a d-2- ry 2- 
symmetry.QB' They may also explain the magneticoLP as 
well as the transporter 1 properties of these systems if the 
gauge fields representing the fluctuations around the MF 
solutions are taken into account. The effect of nonmag- 
netic impurities in the SC state, of. this model has been 
examined by several authors Bcj'tlr 

In this article we study the effect of a magnetic im- 
purity in the superconducting state of the t — J model 
by taking into account both the spin-exchange and the 
spin-orbit interactions between electrons and an impu- 
rity spin. We consider the case of single impurity (lo- 
cated at a site 0), assuming that the concentration of 
impurities is low so that their effect can be treated inde- 
pendently. In this work the impurity spin is treated as 
classical, i.e., the Kondo effect will not be considered. 



The Hamiltonian of our system is 
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where Ci a is an electron operator within the Hilbert space 
excluding double occupancy, and Si and \x are the spin- 
1/2 operator at a site i and the chemical potential, re- 
spectively. For the transfer integrals tij and the antifer- 
romagnetic superexchange interaction Jij , next-nearest- 
neighbor as well as nearest-neighbor terms are taken 



i,j + v 



i,3+T 

with 



t = ±i, ±y and v 



Jij — 



into account 
±i ± y. 

The impurity part consists of three terms: (i) potential 
scattering Vo , which is assumed to be in the unitary limit 
(Vq ^> J,t), (ii) exchange interaction J (> 0) which acts 
between the impurity (Si mp ; Si rnp = 1) and its nearest- 
neighbor sites, and (iii) spin-orbit interaction described 
by H so . The model described by (i) and (ii) is essentially 
the same as the one studied by Poilblanc et alO We 
simplify the model by treating Si mp as a classical spin, 
(Simp)- The spin-orbit interaction H so is given by 
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in the continuum representation, and we transform this 
to that on a lattice by replacing the derivatives with dis- 
crete differences. In the actual calculations we take into 
account interactions only up to second-neighbor sites of 
the impurity for simplicity of numerical calculations. 

We use the slave-boson method to enforce the condi- 
tion of no double occupancy by introducing spinons (fi a ; 
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fermion) and holons (bt; boson) (c; CT = bj/io-), and de- 
couple this Hamiltonian by a mean-field approximation 
(MFA). In the following we consider only the case of zero 
temperature (T = 0), so that holons are Bose condensed. 
Then the mean-field Hamiltonian is written in terms of 
spinons onlyJiJ 
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and the value Vb = 100 is used. For this value of Vo 
the scattering is in the unitary limit, since calculated 
hopping OPs (\) connected to the impurity site vanish, 
indicating that electrons cannot hop to the impurity site. 
The spatial variations of the OPs for Jo = 1 are shown 
in Fig.l. It is seen that is suppressed near the im- 
purity (Fig. 1(a)), and the effect is strongest along the 
diagonals of the square lattice. This is due to the inter- 
ference effects for momenta close to the gap nodes. In 
the region where is not uniform the extended s-wave 
OP (A s ; not shown) and the staggered magnetization 
(Figl(b)) are induced.E^P The important point here is 
that the spin-triplet SCOPs, A px and A py , can also ap- 
pear (Figl(c); Apy has a similar behavior as A px ). A px ( y ^ 

)/2 where 
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where S is the doping rate, and the summations on i, j 
and A are taken over all sites. Here the SC order param- 
eter (OP), Aij, the hopping OP, Xij, an d the magneti- 
zation, (Sf), are defined as 



A« = </i|/ JT >, Xff = (flfj*), 

(s?) = (f^M - /i/a>/2. 



(5) 



The SCOP with the definite symmetry can be con- 
structed from Aij. For example d x 2_ J( 2-wave OP (Ad) is 
defined as 

A d (i) = (A iti+X + A iti - X - A iii+y - Ai,i_ y )/4 (6) 

with A^ = (Ay + Aji)/2 being the singlet-pairing OP. 
Here we do not exclude the spin-triplet pairing, and they 
actually can be finite as will be seen later. 

In a uniform system without impurities only the 
d x 2_ y 2-w&ve OP is finite for the doping rate of our inter- 
est, i.e., 0.15 < S < 0.2 (optimum and overdoped cases). 
In the presence of an impurity, however, OP with other 
symmetries can be mixed and we have to treat their spa- 
tial variations. In order to do this, we diagonalize the 
mean-field Hamiltonian by solving: the following Bogoli- 
ubov de Gennes (BdG) equationEJ) 
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where E n and Ui n are the energy eigenvalue and the 
corresponding eigenfunction, respectively. The unitary 
transformation ^i — Ui n T n diagonalizes the matrix 
Hmfa, and conversely the OPs (Ay, xlj an d (Sf }) can 
be written in terms of E n and Ui n . These constitute the 
self-consistency equations which will be solved numeri- 
cally in the following x3 We use J as a unit of energy 
(i.e., J = 1), and take t/J = 3 throughout in this paper. 

First we examine the effect of spin-exchange scatter- 
ing. The potential scattering is also taken into account, 



(T) 

A\j — (A^ — Aji)/2 is the spin-triplet OP. In the pres- 
ence of the megnetization there is the imbalance of the 
densities of spin-up and spin-down electrons. Then elec- 
tron pairs cannot be formed in singlet channels only, and 
the spin-triplet components appear. We have also exam- 
ined other values of Jo. The results are qualitatively 
the same, and the effect of spin-exchange scattering be- 
comes larger (smaller) with increasing (decreasing) Jo as 
expected. 

The above results can be understood more precisely 
using the Ginzburg-Landau (GL) theory. (The GL the- 
ory is not quantitatively valid except near T c , but it can 
give qualitatively correct results.) The GL free energy in 
the continuum representation can be written asEIP 
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+ K ds [(d x A d )(d x A s )* - (d y A d )(d y A s y + c.c] 
+ K dp [A d {(d x m)A* px - (d y m)A* py }m + c.c] 
+ K sp [A s {(d x m)A* px + (d y m)A* }m + c.c]. 

(8) 

This !Fs is invariant under all symmetry operations for 
the square lattice and we have dropped higher order 
terms. Here we assume that all coefficients in J~s are pos- 
itive except a d . Due to the g d term A^ is suppressed and 
its gradient becomes finite over the range of the coher- 
ence length. Then A s is induced through the mixed gra- 
dient (K ds ) term. The magnetization (denoted as m) is 
assumed to be induced by the spin-exchange interaction 
and have staggered oscillations. When m and Ad coexist 
and have spatial variations, p-wave OPs can occur due 
to K dp term. In Fig.l A px ^ is finite where m is finite, 
and this result is consistent with the above argument. 
All induced OPs are real, since they are determined by 
bilinear coupling terms in Ts- Note the proximity effect 
in a bilayer system composed of a c£-wave superconduc- 
tor and a (anti)ferromagnet can similarly induce p-wave 

opsB 

Next we consider the spin-orbit interaction. It was 
argued, based on a continuum model that the spin- 
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orbit impurity scattering in a d x 2_ y 2-wa,ve superconduc- 
tor can locally create a (d x 2 __aiiiL y )-state which breaks 
time- reversal symmetry (T^ESEiEd) This is because the 
(d x 2_ y 2 + id xy )-state has an orbital moment and its cou- 
pling to the impurity spin can lower the energy of the 
system. We will show that in the t — J model a similar 
T-breaking state can also appear. In order to see this 
we add small t' and J' terms, since the d^-component 
(Aj) is the SCOP defined on the next-nearest-neighbor 
bonds. (Note J' of the order of J is necessary to have 
finite in a system without impurities.) The spin-orbit 
coupling for Ni is estimeted_to be g/a ~30 meV (a be- 
ing the lattice constant). EnEEP Here we use a larger value 
g/a = 2. 4 J in order to compensate the fact that the 
long-range part of H so is neglected in our calculation. 
(J is estimated to be J ~ 0.13 eV, so our parameter is 
about 10 times larger than the realistic value.) We find 
that (Fig2(a), (b)) as well as A s is actually induced, 
but not (S z ) and A p . The induced OPs have nontriv- 
ial phase structure, namely the phases take values other 
than or ir (relative to Ad) and change as functions of 
the position. In particular A^ has a phase tt/2 along 
the diagonal direction where Ad is strongly suppressed. 
This leadSpto a spontaneous current around the impurity 
(Fig.2(c)),E2P and the gap nodes in the d x 2_ y 2-wa,ve SC 
state vanish due to the complex combination of SCOEs. 
These results are consistent with those in ref.l9-21.E-$ 
For the parameters used here the magnitude of the cur- 
rent is of the order of 10 nA near the impurity and is 
smaller in other regions, and the gap due to the induced 
OP is ~ O.lmeV. A rough estimate using Biot-Savart law 
gives a value of ~0.1 G for a magnetic field produced by 
the current at the impurity site. 

Now we examine the case where the spin-exchange and 
the spin-orbit interactions coexist. In this case we find a 
state which carries a spontaneous spin current (Fig. 3) as 
well as usual charge current. The spin current is defined 
as the difference of the currents of spin-up and spin-down 
electrons. Its magnitude is typically 10 -2 of that of the 
charge current. The occurrence of the spin current is 
naturally understood as follows. Spontaneous currents 
can flow because the SCOPs have nontrivial phase struc- 
ture due to the spin-orbit coupling g. In the presence of 
the spin-exchange interaction Jo, the magnetization m 
is finite, indicating that there is the imbalance of the 
densities of spin-up and spin-down electrons. Then the 
contributions to currents from electrons of opposite spins 
are not equal. Hence the spin current arises. 

In summary we have studied the effect of a magnetic 
impurity in the superconducting state of the t — J model. 
The spin-exchange scattering induces spin-triplet SCOP 
near the impurity, while the spin-orbit interaction leads 
to a state with complex OPs and a spontaneous current. 
When both interactions coexist we find a spontaneous 
spin current as well as the usual charge current. In the 
present work the long-range part of the spin-orbit in- 
teraction and the vector potential are not included, the 
latter of which has a feedback effect of the spontaneous 
current on the electronic states.Ei Therefore the results 
presented here can be compared with experiments only 



in a qualitative sense. If some element with a larger spin- 
orbit coupling can be doped in high-T c superconductors, 
there would be a larger (compared with Ni) possibility 
to realize a SC state with a spontaneous current and a 
full gap in the local density of states. The former (lat- 
ter) could be detected by /xSR (low-temperature STM) 
measurements. 
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Fig. 1 Spatial variations of OPs for t = 3, J = 1, 
t' = J' = 0, (i = -0.339 (6 = 0.2), V = 100, J = 1 and 
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g = 0. The system size is 25 x 25 sites, and an impurity 
is located at the center of the system, (a) Ad, (b) A px 
and (c) (S z ). Note that all OPs are non-dimensional. 

Fig. 2 Spatial variations of (a) ReA^, (b) Im A^, and 
(c) the spontaneous current around the impurity. Pa- 
rameters used are t = 3, J = 1, t' = —0.5, J' = 0.2, 
(j, = 0.4 (8 = 0.16), V = 100, J = and g = 2.4. The 
arrows in (c) indicate only the directions of the currents, 
but not the magnitudes. 

Fig. 3 Spin current around the impurity for t = 3, 
J = 1, t' = -0.5, J' = 0.2, (i = -0.5 (6 = 0.189), 
Vq = 100, Jq = 0.6 and g = 4.8. The arrows indicate only 
the directions of the currents, but not the magnitudes. 
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